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Abstract 

There often appear coherently oscillating scalar fields in particle physics moti- 
vated cosmological scenarios, which may have rich phenomenological consequences. 
Scalar fields should somehow interact with background thermal bath in order to 
decay into radiation at an appropriate epoch, but introducing some couplings to 
the scalar field makes the dynamics complicated. We investigate in detail the dy- 
namics of a coherently oscillating scalar field, which has renormalizable couplings to 
another field interacting with thermal background. The scalar field dynamics and 
its resultant abundance are significantly modified by taking account of following ef- 
fects : (1) thermal correction to the effective potential, (2) dissipation effect on the 
scalar field in thermal bath, (3) non-perturbative particle production events and (4) 
formation of non-topological solitons. There appear many time scales depending on 
the scalar mass, amplitude, couplings and the background temperature, which make 
the efficiencies of these effects non-trivial. 



1 Introduction 



Scalar fields play important roles in particle physics and cosmology. The standard model 
Higgs boson, which has been recently discovered, develops a vacuum expectation value 
(VEV) and spontaneously breaks gauge symmetry. In the supersymmetric (SUSY) stan- 
dard model (SM), there are many scalar fields as superpartners of SM quarks and leptons. 
These scalar fields may have large VEVs in the early Universe and their dynamics may 
have significant effects on the baryon asymmetry of the Universe through the Affleck-Dine 
mechanism [T]. 

Inflation in the very early Universe is considered to be caused by a scalar fleld, called 
inflaton, which slowly rolls down the scalar potential [2]. After inflation, the inflaton 
decays into radiation and then hot thermal Universe begins. The inflaton may also be 
responsible for the primordial density perturbation which seeds the rich structure of the 
present Universe. Instead, another scalar fleld, called curvaton, can explain the primordial 
density perturbation [31 lU [5l E] . The curvaton also decays into radiation and fluctuations 
in the curvaton sector are converted to the adiabatic perturbation in the radiation. 

In these scenarios, the understandings of the dynamics of coherently oscillating scalar 
flelds are essential for deriving their phenomenological consequences. For example, the 
inflaton or curvaton must dissipate its energy into the radiation in order for a hot ther- 
mal Universe to begin well before the big-bang nucleosynthesis. It inevitably requires 
interactions between the inflaton/curvaton and other flelds, which somehow will be ther- 
malized. A conventional picture of this reheating process is that the oscillating scalar 
fleld decays into lighter particles and the produced particles are eventually thermalized 
by gauge and/or Yukawa interactions. This simple picture, however, does not necessarily 
hold in general. Depending on the couplings between the scalar (0) and other flelds (x), 
many non-trivial issues appear which makes this topic far more complicated. 

Let us assume the simplest Yukawa coupling between the scalar and a fermion x, as 



in order to induce the (f) decay into radiation, where A is a coupling constant and x is 
assumed to have gauge or Yukawa interactions with SM particles. (As a particular case, 
X itself may be SM flelds.) First, the x flelds obtain the 0-dependent mass through 
this coupling. It is soon recognized that if X(f) ^ m^, where is the mass of </>, the 
decay process — XX is kinematically forbidden except for the region near ~ 0. In 
this case, the perturbative decay rate at the true vacuum does not apply and we must 
take into account particle production phenomena occurring around ~ and subsequent 
thermalization processes of x- Moreover, in general, there always be background radiation, 
or thermal bath, even before the reheating completes. Then, if x is light enough, it is 
thermalized and obtains a thermal mass. If the thermal mass is large enough, the decay 
into X is again kinematically forbidden even if the amplitude is very small. However, 
even if the decay process is kinematically forbidden, can dissipate its energy through 
0-X scattering processes. Thermal bath also modifles the effective potential of through 



^ = A0XX, 
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thermal effects and hence the dynamics of may be significantly affected. 

Therefore, even only the introduction of a simple interaction (11. ip leads to many com- 
plicated issues on the dynamics of (j). Without careful treatment of these issues, we cannot 
discuss cosmological effects of 0. To the best of our knowledge, however, there have not 
been comprehensive and complete analyses on this topic. It was only partly attacked by 
separate literature as follows. 

• Effective potential of the scalar field in thermal environment is well known in the case 
where x takes part in thermal bath [7j. On the other hand, thermal modification on 
the effective potential when x decouples from thermal bath at large value seems 
to be less known except for the context of Affleck-Dine mechanism [S]. However, 
this modification is rather a generic feature if x has a gauge/Yukawa interaction and 
needs to be included in the analyses. 

• Besides thermal modification on the effective potential, also receives thermal dissi- 
pation. Dissipation of oscillating scalar field in thermal bath was studied in Ref . [9] , 
but the applicability of the results are restricted. Ref. [10] considered the dissipation 
of scalar fields in a setup close to ours, although the case of large value was not 
considered. Ref. [H] also studied a similar topic in the case of infiaton. 

• The effect of non-perturbative particle production, or called preheating^ was studied 
in detail in Refs. [12], [131 E] and thereafter. Most of these studies focused on the 
X particle production neglecting the x interaction with thermal bath. Ref. [15] 
considered the case where produced x particles decay into radiation : so-called 
instant preheating. This may efficiently transfers the energy density to radiation 
even in the case where has a large amplitude and the standard calculation of 
perturbative decay into x is not applicable. 

• In association with thermal modification on the effective potential, the scalar field 
may fragment into the non-topological solitons. Once most energy of the oscillating 
scalar field is absorbed by solitons, estimation of the dissipation rate and the decay 
temperature of may be significantly modified. This is studied in the context of Q- 
balls when the complex scalar oscillates with an elliptical orbit in the complex plane. 
The formation of solitons, however, can occur even in the case of a real scalar. Its 
phenomenological consequences have not been investigated in detail in a setup we 
are interested in. 

In this paper, we consider the dynamics of having interaction of the type (11. ip in 
thermal environment. We take into account various effects listed above in reliable manners 
and study in which situation these effects become substantial for the scalar dynamics. As 
for the issue of non-topological solitons, we avoid to give definite conclusions but mention 
possible consequences. We only have four model parameters : m0,0j,A,a(= g'^/Air), 
where is the bare mass of 0, 0j is the initial amplitude of 0, A is the Yukawa coupling 
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between and x, g is the coupling constant between x and thermal bath, which is typically 
an SM gauge coupling. In addition, we introduce the reheating temperature Tr, which 



controls the density of thermal bath|^ Even in this simple setup, there appear many 
time scales and we must be cautious about the usage of a particular formalism case-by- 
case. Eventually, we have successfully figured out the global picture of the scalar dynamics 
in broad parameter spaces. We believe that, although this is a simplified toy model, it 
captures essential features of the most applications to realistic models. This is because 
oscillating scalar fields must decay into radiation unless it becomes (a part of) dark matter, 
and then it inevitably couples to lighter species, which in turn have interactions with SM 
particles, and this is nothing but a situation we are considering. Even if couples to 
many fields with coupling strength varying by orders of magnitude, the dynamics mainly 
depends on the largest coupling and hence our analyses apply. 

This paper is organized as follows. In Sec. [2] we describe our basic setup and list some 
particle physics motivated examples. In Sec. Owe review thermal effects on the scalar field 
evolution based on the closed time path formalism. In particular, thermal modification 
on the effective potential and the dissipation rate will be introduced, which significantly 
affects the scalar dynamics. In Sec. HI the effects of non-perturbative particle production 
are discussed. This also serves as a dissipation of the oscillation of the scalar field. In 
Sec. [SI we solve the scalar dynamics taking these effects into account. Readers who are 
only interested in results may skip preceding sections and only read Sec. El We conclude 
in Sec. El 



2 Overview 
2.1 Setup 

Let us start with a model where a real scalar couples to another field x (fermion) and 
X (boson) at the renormalizable level. The Lagrangian is given by 

= C^in - \ml<P^ - A(j)xx - \h^^x^ - X<PxX: (2.1) 

where A, h and A are coupling constants and £kin denotes canonical kinetic terms. The bare 
mass for xix) is neglected in what follows. An additional assumption is that the xix) has 
gauge or Yukawa couplings, which is collectively denoted by g, to particles in thermal bath. 
The introduction of xix) field in (12. ip should be regarded as schematic. Our arguments 
are not affected even if there are many fields which couples to and the interaction term 
is given by £ = ^ij^PXiXjy etc. For example, if xix) ha,s gauge interactions, summation 
over the gauge indices should be promised. Notice that in a non-SUSY theory, there is a 

In this paper we assume the existence of thermal bath, which comes from the inflaton decay. Thus 
our analysis of the scalar dynamics is not applied to the inflaton itself. The inflaton dynamics is also an 
interesting subject, which will be discussed elsewhere. 
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large correction to the scalar potential. We have the Coleman- Weinberg potential [16] as 



V- m^(0) r m2(0) 3 

F LA*. 

where ei;' = +1 for a real scalar boson x, and ep = —2 for a (Weyl) fermion x- This yields 
the quartic potential of the form Vcw ~ In 0, In 0. It must be included in the 
analysis of scalar dynamics unless the couplings h and A are small enough to be neglected. 

In SUSY, and x are components of chiral supermultiplets. In this case, assuming 
the superpotential of the form 

W = A0XX, (2.3) 
where we used same symbols for the chiral superfields, we obtain the Lagrangian 

£ = £kin - (A0XX + h.c.) - {\X? + - Ax?\t? - ^SB, (2.4) 

where £kin denotes the kinetic terms of the fields, x{x) ^i-nd xix) ^'^^ ^ fermionic and scalar 
components of each chiral superfield. In SUSY, therefore, the Yukawa coupling and the 
four-point coupling are both given by A. In all the analyses, we assume this for simplicity. 
Note that scalars are complex fields. SUSY breaking effects induce scalar potential for 

as 

VsB = m2|0|2 + mJIxP + m||^|2 + (A^0xx + h.c.) , (2.5) 

where m^,m^,m^, are SUSY breaking parameters of same orders of magnitude. One 
of the good points for considering SUSY is that radiative corrections to the scalar po- 
tential gets suppressed and we do not need to worry much about the Coleman- Weinberg 
correction, since quartic terms in cancel out and there only remain small logarithmic 
correction to the quadratic terms from the SUSY breaking effect. We implicitly assume 
this in the following. 

After all, in order to avoid unnecessary complexity, we consider a following phenomeno- 
logical model 

C = ^kin - ^m,|02 - A0XX - AV'x', (2.6) 

where is a real scalar, x is a fermion and x is a scalar boson. We do not include scalar 
tri-linear interactions because its typical size is much smaller than others. We study the 
dynamics of scalar field in the presence of thermal background interacting with x and x. 
We sill see that the resulting scalar dynamics in this simple theory is highly non-trivial. 
Hereafter, the coupling constant A and the x coupling with thermal bath, y/Ana), is 
assumed to be smaller than unityl^ 

Let us suppose that has an initial value of 0j during inflation and follow the dynamics 
of 0. Conventional arguments are as follows. After inflation, is frozen at 0j until the 
Hubble parameter decreases to m^, at which begins to oscillate with an initial amplitude 

For notational simplicity, we often call "x" both for fermion x ^-nd boson x unless we need to 
distinguish them. 



(2.2) 
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of 0j. The abundance of 0, in terms of the energy density (p</,) to entropy density (s) 
ratio, is then given by 



where Mp is the reduced Planck scale, Tr is the reheating temperature after inflation and 
Tos = {n'^g:^/90)~^^^^yn^Mp. Then, decays into x with the decay rate given by 



r0(0 xx) 



3 /2 

^fl-M^I . (2.8) 



i7i \ 



Since x is assumed to have sizable couplings to thermal bath, the energy stored in the 
form of coherent oscillation goes to thermal bath after the decay. 

The above arguments, however, miss some possibly important effects. First, couples 
to thermal bath through x particles. Although x itself may be heavy enough to decouple 
from thermal bath at large value, it may still modify the effective potential of 0. Second, 
if the finite-temperature correction to the scalar potential is dominant, the fluctuation 
develops to form non-topological solitons. Once fragments into solitons, their cosmolog- 
ical evolutions are significantly different from the case of homogeneous oscillation. Third, 
the perturbative decay rate (12. 8p is not simply applied since x can be heavier than 
during the course of oscillations with large amplitude. Instead, non-perturbative particle 
productions when passes through the origin may be efficient, which in turn reduces the 
energy of coherent oscillation. 

Therefore, even in the simple model like (12. 6p . the dynamics of scalar field is compli- 
cated. There appear many time scales related to the above mentioned effects, depending 
on the parameters m^, 0j, A and the background temperature T. Without taking into 
account these effects, one cannot evaluate the the abundance of and its cosmological 
consequences. In the following sections, we discuss the dynamics of in detail. 

It should also be noticed that the following analyses are applied to the more general 
forms of the scalar potential for 0. For example, let us assume that the scalar potential 
takes the form 

V{<f)) = -ml<p^ - -H'<P^ + (2.9) 

where n > 4 is an integer, H is the Hubble parameter and c is an 0(1) positive constant. 
In this case, the field sits at = {cH'^Mp~'^/ kY^^"'~'^^ for H > m^, and then it begins to 
oscillate around the quadratic potential with initial amplitude of 0j = {cm'^Mp''^/ k)^/("~^) 
at if ~ if there are no thermal corrections. It is not difficult to apply our results to 
such a case. 
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2.2 Examples 

Before going to discuss details of the scalar dynamics, we point out that broad class of 
models discussed in particle physics and cosmological contexts are actually described by 
(EEl). 

2.2.1 Curvaton 

Curvaton is a hypothetical scalar field, which is responsible for the generation of observed 
cosmological density fluctuations [31 IH |5l [6] . In the simplest case, the curvaton (0) has 
a quadratic potential and is supposed to be frozen at the initial value (j) = (pi during 
inflation. Since the curvaton mass is mush smaller than the Hubble scale during 
inflation, -ffinf, it obtains quantum fluctuations characterized by -ffinf, which in turn seeds 
the density perturbation with nearly scale-invariant power spectrum. Since the curvaton 
must decay in order to convert the fluctuation in the curvaton into the radiation, we 
need to introduce couplings between curvaton and some particles, x? which interacts with 
standard model particles. Thus the model of (12. 6p . with/without 0^%^ term, describes the 
'minimal' curvaton model. 

An interesting feature of the curvaton model is that it can generate a large (local- 
type) non-Gaussianity in the cosmological density perturbation, while it is difficult in the 
standard single-field inflation model [171 [HI [iSj- The curvature perturbation ( in the 
curvaton model is given by [H] 

C = Cg + ^/NLC|- (2.10) 
Here Q denotes the Gaussian part of ( and is given by Q = -RiJinf/(37r0.j) where 

3P0 



R 



3p^ + Apr 



(2.11) 

ip decay 



which roughly describes the energy fraction of (p in fhe Universe at the decay. The non- 
linearity parameter, /nl, is evaluated by 

. 4 (l « _ (2,2) 



4R\ 3 

It is clear from this expression that we need R < 0.1 for obtaining /nl ~ (9(10). This 
means that the curvaton must not dominate the Universe before it decays in order to obtain 
large non-Gaussianity. Therefore, the initial condition should satisfy 0j -C Mp, since 
otherwise the curvaton eventually dominates the Universe unless the inflaton reheating is 
so late. However, as will be seen later, thermal effects are always important if (pi <C Mp. 
Then the estimates of the curvaton abundance and hence the resulting non-Gaussianity 



are significantly affectedPn Thus it is important for the curvaton scenario to reconsider 
the scalar dynamics. 



The estimate of the non-Gaussianity is affected due not only to the change in the curvaton abundance 
but also to the departure of the scalar potential from the quadratic one [2^1 EH • 
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2.2.2 Flat directions in MSSM 



Next, let us consider a supersymmetric (SUSY) theory. In SUSY, there appear many 
scalar fields as superpartners of the standard model fermions. It is known that, in the frame 
work of the minimal SUSY standard model (MSSM), some combinations of sfermions have 
vanishing F- and D-term potential at the renormalizable level in the SUSY limit [H[22|[23]. 
Such a flat direction in the MSSM, parametrized by a complex scalar $, obtains scalar 
potential from the SUSY breaking effect as V{<^) = m|,|$p and also from possible non- 
renormalizable terms in the superpotential or Kahler potential. The dynamics of flat 
directions was investigated in detail in the context of Affleck-Dine baryogenesis [H [22] 
including thermal effects [Ml El |25] and the formation of Q-balls [261 [271 il |29l [301 EH ESI 
[33] . While the motion of flat direction in the complex plane must be elliptical in order to 
generate sizable baryon number, it is possible that the motion is nearly one-dimensional if 
the y4-term contribution to the scalar potential is somehow small, as in the gauge-mediated 
SUSY breaking models. Then the dynamics is essentially one dimensional described by 
the motion of the radial component of $. Since flat directions, which are composed of 
squarks, slepton and Higgs bosons in the MSSM, have gauge and Yukawa interactions, 
the interaction terms like (12. 6p appear. In this case, x's are MSSM flelds and A's are 
corresponding gauge or Yukawa couplings. Thus the dynamics of the flat direction also 
reduces to a simple model of (12. 6p . 



2.2.3 Right-handed sneutrino 

Right-handed neutrinos are often introduced because it can explain the tiny left-handed 
neutrino masses through the seesaw mechanism, and also the observed baryon asymmetry 
of the universe through the leptogenesis scenario by their non-equilibrium decay. In SUSY, 
there are scalar partners of them, right-handed sneutrinos denoted by Ni {i = 1, 2, 3). The 
superpotential is given by 

W = im^^iV.iV, + yf^N.LjHu, (2.13) 

where Ni, Lj and Hu denote the right-handed neutrino, left-handed lepton and up-type 
Higgs superflelds, mjv, is the Majorana masses for the right-handed neutrino and is the 
Yukawa coupling among them. Focusing on the lightest component, A'^i, and writing the 
D-flat direction of the LiHu as LiHu = x^, the scalar potential is given by 

V = \m^,N, + yfixT + VsB. (2.14) 

Here Vsb denotes the SUSY breaking effect including terms such as BmN-^NiNi+h.c. with a 
constant B of the soft SUSY breaking scale. Assuming that initially A'^i has large amplitude 



and X = Oj^the scalar potential of the radial component of Ni is simply quadratic and, as 



a result, the dynamics is described by (12. 6p . Cosmological implications of the right-handed 



This assumption is valid if mjvi ^ -ffinf ■ 
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sneutrino condensation, mostly in the context of non-thermal leptogenesis, were discussed 
in the literature [SU [351 IMl ETJ |38l [391 HO]- Thermal correction on the right-handed 
sneutrino dynamics was partly discussed in Refs. [371 1321, but the complete analyses on 
this system were not performed so far. It may be crucial for the estimate of the baryon 
number generated by the right-handed sneutrino decay. 

2.2.4 Peccei-Quinn scalar 

The most attractive solution to the strong CP problem is the Peccei-Quinn (PQ) mecha- 
nism [in 1^ . In the so-called hadronic axion model [13], the complex PQ scalar field, $, 
couples to vector-like quarks Q and Q as C = X^QQ. If the PQ scalar obtains a vacuum 
expectation value of order of 10^-10^^ GeV, the strong CP problem is solved without con- 
flict with cosmological/astrophysical observations related with the axion phenomenology. 

In SUSY, there is a fiat direction in the PQ scalar sector, called saxion. The saxion 
feels the scalar potential from the SUSY breaking effect. Thus the saxion potential at 
large field value is often quadratic, although this is not mandatory. The dynamics of 
the saxion and its cosmological implications were discussed in Refs. [lH [lH [l6l [13 [HI 
[13 [Sni [SB [S21 [S31 [51 [23 [SSI [SZ]. if the initial amplitude of the saxion is much larger 
than the PQ scale, the saxion dynamics resembles the model described by fl2.6p . where x 
identified with the PQ quarks, Q and Q. (Thermal correction on the saxion potential in 
this context was pointed out in Ref. [SH [SZ]-) One should notice that the final stage of 
the PQ scalar dynamics toward the symmetry breaking minimum may include additional 
steps. For example, if the saxion is trapped at the origin due to the particle production 
effect, the saxion finally causes thermal inflation [5H1 [59]. However, the first stage of the 
saxion dynamics still falls into the model of (12. 6p . Details of the saxion dynamics certainly 
depend on how the PQ scalar is stabilized, and more complete analysis will be presented 
elsewhere. 

3 Thermal effects on the scalar field dynamics 

In this section we discuss thermal effects on the dynamics of homogeneous scalar conden- 
sation in this model. Since the field interacts with thermal plasma through x fields, 0's 
dynamics may be affected. In particular, there are two effects caused by thermal plasma: 
thermally modified potential for and the dissipation of </> to thermal plasma. Thermal 
effects are qualitatively different whether the field amplitude is large or small: \4'{t) ^ T 
or \(j){t) <C T. Hence, we study these regimes separately in the following. Basic ingredients 
for closed time path formalism are summarized in Appendix [C] briefly. The expectation 
value of operator A with respect to the physical state p is denoted by 

(i) := tr{pi}. (3.1) 

In the following, we assume that the system has a spacial translational invariance and the 
homogeneous condensation of ip is denoted as = (ip) . 
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3.1 Large Amplitude Regime 

First, let us study the large amplitude regime: X(f){t) ^ T. We separate the relevant time 
scales and follow a coarse-grained effective motion equation for </>. 

As shown a posteriori, a typical dynamical scale of scalar condensation (f) at this regime 
is much slower than the thermalization time scale: <C Tth, if the cj) oscillates with the 
thermally modified potential. Of course, this condition is automatically satisfied if the 
oscillates with the vacuum mass m,^ with <^ Fth- Then we can assume that the 
"fast" fields in thermal bath feel the "slow" 0's dynamics as almost static. Hence, on the 
typical time scale of 0's dynamics, the other fields in thermal bath have long enough time 
to thermalize with an each value of background field and completely forget about their 
past. In such a regime, the equation of motion can be simplified by tracing out all the 
degrees of freedom in thermal equilibrium. 

There is one more important aspect at this regime. Since the field amplitude is large 
X(f){t) ^ T, the number densities of x particles are very small. Therefore (j) does not feel 
the existence of these particles directly. However, this does not immediately mean that the 
(/)'s dynamics is completely free from thermal plasma, despite the fact that the interaction 
between (j) and thermal bath is only mediated by the heavy x fields. This is because the 
free energy of thermal plasma depends on the background (p field at a higher loop order. 



and hence the "pressure"!^ dominantly affects the 0's dynamics if the vacuum mass 
is very small. 

With the above consideration in mind, let us derive the effective equation of motion of 
0's homogeneous condensation. Initially, the system is prepared as the canonical ensemble 
with the non- vanishing background field expectation value (j){ti). In principle, one may 
follow all the dynamics in terms of equations of motion derived from Closed Time Path 
2PI (nPI) effective action: = 6T[(f),Gij]/6(f), = 6T[((),Gij]/SGij where the subscript i 
runs all the fields, but it is practically difficult (See Refs. [601 EI] and references therein). 
Since all the fields except remain in thermal equilibrium with the background field 
we can safely assume that the propagators Gij are given by thermal ones. Therefore, the 
effective equation of motion is reduced to 

6(f) 6(j) 6(j) 

where is the 0's classical action and F denotes the sum of bubbles calculated in terms 
of the thermal propagators with the background field 0. 

Let us evaluate F approximately. Since the motion of is very slow compared to 
the thermal plasma, we will neglect the time dependence of as a first step, and then 
take it into account approximately. If the is regarded as static, F is merely the "free 
energy" of thermal plasma with the constant background field 0, and hence we have 
(5F/(50 ~ —dVes/d(j). Since the x particles are absent at T ^ X(f){t) due to the Boltzmann 



Analogous to P = — ^ in statistical mechanics. 

ox 
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suppression, we can integrate out x fields first in calculation of the free energy. This leads 
to the effective operator which contributes to the running gauge coupling constant g 



A 



ln(A2(/)7T2) f^^'^F;^ (3.3) 



where A is a constant determined by the representation of x fields. Because the free 
energy of hot QCD plasma has a contribution proportional to g^{T) T^, this term induces 
the so-called thermal logarithmic potential [S]: 

Veff D aa{Tf ln(AVV^') (3-4) 



#6 



where a is an order one constant 

The typical dynamical scale of this term is given by JTi^g^ ~ a^T^/cf)^. Thus, if the (p 
oscillates with the thermal logarithmic potential, the motion of </> is adiabatic with respect 
to the typical thermalization time scale of thermal plasma, Fth ~ aT ^ "i^g, at the large 
amplitude regime, T ^ X(f). 

Next, let us take into account the time dependence of (j) at the leading order, and derive 
the dissipative coefficient, which describes the typical relaxation time scale of (p. We will 
separate the field into the constant background v at time t and the small deviation 
and then evaluate T up to the first order in Scpit). T can be expanded as 



srn+li 



1 Sf[v + S(p] 

V 6(j){t) V^nlJ^"" 6(f){t)6(f){ti)---6(j){tn) 



5(P{h)---5(P{Q (3.5) 

(j)=V 



d(f) 



y"dr nret(t-r,0) 50(r) + --- (3.6) 



up to the first order in S(j). Here V denotes the spacial volume. The first term is the free 
energy, which is already obtained above. At the leading order, the self energy is given by 

n,et(x) := -^^(x°)nj(x); Uj{x) = ([6(x), 6(0)]) (3.7) 

where the effective interaction, obtained from integrating out x fields, is given by 

Ant = 500; 6 := -^F^a^'^F;,. (3.8) 

Here the ensemble average (■ ■ ■) is merely the thermal one, because the thermal plasma 
remain in thermal equilibrium. Since the relevant time scale of the self energy is determined 



The sign of the coefficient a depends on the modcL In the model we arc considering, % is a matter 
field, which is either a fermion or a scalar boson having gauge interactions. In this case, a is positive. On 
the other hand, if % is a gauge boson, a can be negative. In all the analyses in this paper, we assume that 
o is positive. 
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by thermal degrees of freedom, it is much faster than one of 0's dynamics. Then we can 
approximate 50(r) as 50(t) ~ 4>{t){T — t), and the following equation is obtained 

dr n,et(t - T, 0) 5(P{r) = -i j dr 53n,et(r, 0) r 

= -lim ^"-^(^'°^ m (3.9) 

In the first equality, we have used the fact that the real and imaginary part of Ilrctit^ 0) 
are even and odd functions in t respectively. As can be seen from Eq. (13. 9p . the dissipative 
coefficient is imprinted in the imaginary part of self energy: 

r, := - Ita = ,™ !k(^. (3.10) 

In the second equality, we have used the Kramers- Kronig relation. 

Note that the operator O is related to the trace anomaly of gauge field. As pointed 
out in Refs. [631 IM] , the dissipative coefficient induced by the trace anomaly is directly 
related to the bulk viscosity of hot QCD plasma [65] : 

C = ^ Imi ^/ d'x e"-*([T%(t,x),T^(0,O)]), (3.11) 



where 



and it is estimated as 



T%(^,x)^-|f"'^'^F- &o:=^|^, (3.12) 



A V 36C (Ay (127ra)2 
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Here the bulk viscosity C, is evaluated at the weak coupling regime 

Finally, with taking account of the adiabatic expansion of the universe H <^ Fth, we 
obtain the effective equation of motion of at the large amplitude regime: 

0+(3ff + r^)0 + m20 + V^g = O (3.14) 

where the dissipative coefficient and the effective potential are given by 

[s^-J lna(T)-i ^ ^^-^^^ 
Veff(0)~ aa(T)2TMn(AV/T') (3.16) 

respectively. Here we have omitted the constant term of effective potential which is inde- 
pendent of the (j) at the large amplitude regime. 



T^'' Here we assumed iriq ^ aT where nig is the heaviest zero temperature quark mass in thermal bath. 
"In thermal bath" means that the zero temperature mass is at most T. 
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3.2 Small Amplitude Regime 

Second, let us study the small amplitude regime: X(j){t) <^ T. In this regime, x's number 
density can not be neglected and the x particles in thermal plasma may directly affect 
the 0's dynamics. In the following, we assume that the dynamics of </> is not so violent as 
x's number density cannot remain the Bose-Einstein distribution. In other words, the x's 
propagators can be well approximated with thermal propagators. Typically, this is the 
case where 0's dynamics is slow enough for thermal plasma to remain thermal equilibrium. 
A possible non-perturbative production when the passes through the origin is discussed 
in the next Sec. H] and the applicability of discussion given in this section is clarified. 
In the small amplitude regime, the effective equation of motion can be expressed as 

+ 3H<P{t) + mj0(t) + J dr U,,,it - r, 0) 0(r) = 0. (3.17) 

The self energy can be decomposed into the local and non-local parts as 

nret(t, 0) = M^5{t) + nret(t, 0), (3.18) 

where 

= 2\\i'), (3.19) 

nret(t,o) = -ie{t)iij{t,Q)- nj(t,x) = ([6(t,x),6(o,o)]), (3.20) 

at the leading order in A and A, with the operator O being 

d = \xx + ^XX- (3.21) 



Here the ensemble average (■ ■ ■) is also the thermal one 



#8 



The self energy is computed perturbatively in terms of thermal propagators of x 



fields!^ First, let us evaluate the local part and the dominant real part of self energy. 



which lead to the thermal mass proportional to T^. The local part is given by [75] 

= / - — - ■'^^ ~ , 3.22 

J 27r^ u;k 6 ' ^ ^ 



where = \/ + k^, at the leading order in X(j)/T, with fs being the Bose-Einstein 
distribution. Similarly, the dominant real part of self energy, which is caused by the 
Yukawa interaction, can be evaluated as 

afjfiret ^ (3.23) 



Although we do not consider scalar tri-linear coupling in the application in the following sections, 
here we give formulae for the case of scalar tri-linear coupling just for illustration. 
See Appendix. IC.2l for explicit forms of thermal propagators. 
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at the leading order in X(j)/T [66]. Note that if the amphtude is large, then these "thermal 
masses" of are absent due to the Boltzmann suppression of x particles. Hereafter thermal 
masses are denoted as m^j^(T) collectively where the superscript i denotes the species: (j), x- 
Next, let us evaluate the dissipative coefficient imprinted in the imaginary part of self 
energy. The dissipative coefficient is given by 



Sriretlw, 



nj(w,o) 



2uj 



(3.24) 



where the effective mass of 0, m'^g, is given by 



m 



eff 



m. 



(3.25) 



at this regime. If the effective mass is small compared to the typical thermalization time 
scale, the dissipative coefficient can be well approximated with — )■ 0, as in the case of 
the oscillation with the thermal logarithmic potential. The self energy can be expressed 
as 

(/^(go) - fpiqo + "^eff)) tr pJ(go, q)pj("^fff + go, q) (3.26) 



+ A2 



(/i?(go) - fBiqo - mfg) j pJ(go, q)pf (mfg - go, q), (3.27) 



where is the spectral density for the bosonic/fermionic x ^^id and fs/p is the Bose- 
Einstein/Fermi-Dirac distribution respectively. 

For simplicity, we will assume that the spectral density is well approximated by the 
Breit-Wigner form. The Breit-Wigner form of spectral density for boson is given by 



Py(go,q) 



2gorc 



[q'o - ^l? + [gorj2 



(3.28) 



where Qr 



m^^{TY + A^0^ + q^. And the Breit-Wigner form for fermion is given 



by W. 

pJ(go,q) = 



s=± 



[go - %f + r^V4 



(7o - q ■ 7) + 



q 



[go + %Y + r^74 



(7o + q ■ 7) 



(3.29) 



Here the plus (s = +) and minus (s = — ) contribution correspond to the ordinary particle 
like excitation in vacuum and the new collective excitation in thermal plasma, so called 
the plasmino, respectively [691 EH], is the quasi-particle energy and Fp corresponds 
to the typical relaxation time scale of quasi-particle, so called the thermal width. In 
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the following, for simplicity, we will neglect the plasmino contribution and the dispersion 

relation is approximated by fip ^ \J m^^{TY + where ra^ is the asymptotic mass, 

which is given by ra^ = \/2m^{^ [68]. Note that Eq. fl3.29p is valid at the very small 0's 
amplitude regime. If the 0's amplitude is not so small T ^ A0 ^ gT, then the spectral 
density of fermion can be approximated with [10] 



pJ(go, q) = (A0 + i) [,g_4p + [gor,]^ ' ^'-'^^ 

where cUq = \J X^(j)^ + q^. 

To evaluate Eqs. fl3.27p and f l3.26p . let us consider typical two cases: (i) the 0's dy- 
namics can be regarded as adiabatic with respect to thermal plasr na, and (ii) the 0's 



amplitude is small Xcf) <^ gT where (p is the maximal amplitude of l''^^"! (not necessarily 
In the first case (i), m^g ^ aT, the dissipative coefficient can be approximated with [10] 

1 2 /• j4 



2T 



J (2^'^^*^^°'* ~ [px(^o, q)pJ(go, q)] (3.31) 

A2 r d^q 
2fJ (2^ 



+ 7^ / TZ^fBiQo) (1 + fsigo)) P?(go, q)p?(go, q) (3.32) 



Using the Breit-Wigner form for the spectral density, one finds the dissipative coefficient 
from Eq. fl3:32ll as 

r? ~ ^, (3.33) 

where the thermal width is roughly approximated by Fq ~ aT. On the other hand, the 
dissipative coefficient from Eq. fl3.3ip can be roughly evaluated as 



X^aT for A0 < aT 

X^-^ for aT<A0«T, 
al 

where the thermal width is roughly approximated by Fq ~ aT. Note that X'^x^(l)'^ term 
also leads to the same order contribution as the latter one. In addition, note that the 
result in the interval between two regimes, A0 ~ aT, is a rough approximation, and 
that the latter expression is not applicable at A0 ~ T since the exponential suppression 
factor in the integrand dominates. As pointed out in Ref. [71], above one-loop results are 
merely approximate ones because higher-loop contributions such as ladder diagrams are 
comparable to one-loop ones with the vanishing external energy w — t- 0. This is much like 



"^^^ In this case, the efficient non-pcrturbative particle production is absent. See also Sec. |4] 
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what happens in the calculation of viscosity coefficients from Kubo formulas [TQl [65]. As 
discussed in Ref. [10], the resummation of infinitely many diagrams can change the one- 
loop result by several factors. In the following, however, we will estimate the dissipative 
coefficient with one-loop diagrams as a rough approximation. 

In this case (i), the obtained equation is formally equivalent to Eq. (13.141) . since the 
free energy has the thermal mass term X^T'^cfP' at the small amplitude regime. 

In the second case (ii), the dissipative coefficient is easily obtained if the "decay" of (f) 
to X is kinematically allowed |llj : 



r? ^ T^\h- ^ "^K^Y + 2/b«/2)) Oimt, - 2m^^{T)) (3.35) 



" ^ V ' - '^^f^ - 2/HrnJff/2)) e{mt, - 2m-^(r)). (3.36) 



If the decay is kinematically forbidden [aT <^) m^g ^ gT, then the dissipative coefficient 
is from the tail of Breit-Wigner distribution due to the thermal width, and hence it is 
suppressed by the coupling in thermal plasma g, compared to the "decay" [72| [TT]. 



4 Non-perturbative particle production 

Let us discuss the effects of non-perturbative particle production in this model. Since the 
scalar field (f) passes through the origin during the coherent oscillation, the adiabaticity of 
the coupled particles x is necessarily broken. It implies that the x particles are produced 
in each oscillation even if the perturbative decay of cf) is not efficient. This is called the 
preheating [T2| [T^. For a boson x, the production efficiency is enhanced as the x particle 
number increases and this leads due to the parametric resonance effect if the dissipative 
effect of X is not large [77]. Even for a fermionic Xi the preheating can have a significant 
effect on the reduction of the energy density of cf) through the instant preheating [TJ], if 
X has sizable interaction with other particles. (See Refs. [75| [7^ for a theory of fermionic 
preheating.) 

4.1 The case of zero temperature mass 

First, we consider the case where the cf) oscillates with the zero-temperature mass and 
its time dependence is expressed as (/)(t) = (psmim^t) in one oscillation. Most discussion 
below applies for both fermion x and boson x as long as their coupling constants to are 
same (see (12.61) ). and we do not distinguish them unless otherwise stated. The mass of x 
varies with time and it has a frequency 

ujl = e + mUTf + \^{t)\ (4.1) 
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where /c is a wavenumber and m^^{T) denotes the thermal mass of x? given by 

mliT) ~ gT, (4.2) 

Notice that it is not evident that we can take rri^^ ~ gT since the (f) oscillation frequency 
may be larger than the x thermalization rate. However, we will see that in the practical 
use, Eq. (14. 2 p gives appropriate results. 

Particle production occurs if the following adiabaticity condition, |a;^/a;^| <C 1, is 
violated [H], where we have 



a;2 

X 



A^0^m<^ sin(m^t) cosim^t) 



[F + m^h(T)2 + sin2(m<^t)]3/2 ' 
First, consider the case of m^^{T) ^ m^. The typical wave number k^, below which the x 



mode is amplified, is given hj = y Xm^cp. The typical time interval At^, in which the 
adiabaticity is temporary violated, is estimated to be At* ~ This is much smaller 

than the oscillation period for <^ A0. The created x particle number density during 
the passage of the minimum ~ is given by 



^^^(A™^ ,44) 



Next, consider the opposite case : m^j^(T) ^ m^. The condition |cj^/a;^| ^ 1 requires 
kl ^ m^j^(T)2, which is rewritten as Xcp ^ m^^{Ty/m^. As long as this condition is 
satisfied, a typical time interval, in which the particle production occurs, is much smaller 
than the oscillation period : At^, <^ l/m^. Thus we impose the following condition for 
efficient particle production : 

r g-^T^ ] 

\(f) ^ max < m^, > . (4.5) 

I J 

If this condition is met, modes with A; < /c* is amplified and the produced number density is 
estimated by Eq. (14. 4p . Note that this implies that if the amplitude is very small <^ gT, 
the efficient non-perturbative particle production does not occur. 

Subsequent evolution of the system crucially depends on the decay /dissipation rate of 
X • ^3- After the passage of ~ 0, the x mass increases and correspondingly the decay 
rate of x also becomes large. Assuming the typical x decay rate as ~ am^ = aX\(f)(t)\, 
it decays at tdec ~ {o!Xm^(j))~^^'^, well before the again reaches the maximum at the 
opposite side of the potential, if ^ aX(f). Otherwise, returns back to = before x 
decays into radiation. Then the parametric resonant amplification of the x modes occurs 
for a bosonic x- For a fermionic x, the Pauli blocking suppresses the further particle 
production. Note that x also has a thermal dissipation rate of order of F^ ~ aT for 
A0 < T, but this does not dominate the above estimate as long as the condition (14. 5 p is 
met. 
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(a) In the case of A0 ^ m^/a, a phenomenon similar to the so-called instant preheat- 
ing [15] takes place. In each oscillation, the x particle number density ( 14. 4 p is produced, 
which soon decays into radiation. Therefore, the fractional energy density which loses 
in one oscillation is given by 

^ ^ (4.6) 

As long as the preheating lasts, this fractional amount of energy dissipates into thermal 
bath in one oscillation. The resultant energy density is given by 

pS) - P*.(l - ^^T'"'' = - hT'"^- (4.7) 

Unless is exceedingly small, the energy density efficiently is transformed into the 
radiation and the amplitude decreases correspondingly. The effective dissipation rate of 
is then given by 



„^„„<p- ^ ^ ^. (4.8) 

Notice that this dissipation exists in the regime the perturbative decay of into x is 
prohibited, i.e., when viv^i^ ~ A^0^ + g^T'^ ^ m^. This regime ends when the amplitude 
decreases and either of the following condition is satisfied, (i) A0 = 47raT^/m0, where 
thermal mass of x becomes efficient. Then the particle production events end. (ii) A0 = 
m^/a, where the regime of instant preheating ends. Then it enters the regime (b) in the 
following. 

(b) In the case of A0 <^ m^/a, the produced particles around ~ 0, which fills the 
phase space density of /c < /c*, survive until the again returns back to ~ 0. For a 
boson Xi the efficiency of particle production increases due to the parametric resonance 
effect. As a result, the x number density exponentially increases until the backreaction 
terminates the resonance. It is expected that at this stage the and x energy density are 
equilibrated. For a fermion x, on the other hand, no such resonant effect takes place due 
to the Pauli blocking effect. In either case, at most 0{1) fraction of the energy density 
is dissipated. 



4.2 The case of finite temperature mass 

Next, let us consider the case where the oscillates by the finite-temperature effect, and 
see whether or not the non-perturbative particle production occurs when the passes 
through the origin. 

First note that the rate of thermalization in thermal bath, Fth ~ aT, is larger than the 
effective mass scale of at the beginning of oscillation, Fth ^ ~ ^ts- This observation 
immediately means that if the begins to oscillate with the thermal mass potential, then 
the motion of is adiabatic with respect to the thermalization time scale of thermal 
plasma, mfj^(T) ~ AT <^ Fth ~ oiT, as far as it oscillates with the thermal mass. Hence, 
the discussion given in Sec. 13.21 is applicable at all times as shown below the case (i). 
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Second, if the oscillates with the thermal logarithmic potential (not necessarily at 
the beginning of oscillation), then the typical thermalization time scale of thermal plasma 
is larger than the effective mass scale of </>, mff^ ~ aT^/cj), because of > (^c = T/X. 
Thus the discussion given in Sec. 13.11 is applicable in this case, and we can reliably use 
the thermal logarithmic potential for (pc < <P{i)- At (f){t) < (pc the x particles begin to be 
produced from thermal plasma and they yield a finite density correction to 0's equation 
of motion. 

At this stage, we have to distinguish two cases, 
(i) A < a 

In this case, the motion of is adiabatic with respect to the relaxation time scale of 
thermal plasma. Thus, the discussion given in Sec. 13.21 is applicable. As discussed in Sec. 
13. 2[ the finite density correction can be expressed as 

J dr n,et(t - r)0(r) ~ mfjT) ^ + (4.9) 

where the thermal mass is given by ra^^{T) ~ AT. Then, let us see whether the non- 
perturbative particle production occurs or not when The effective mass of x is given by 
m^g^ ~ aT^ + A^0(t)^. The adiabaticity parameter is estimated to be 



X 



A^c A 

< ^ = - « 1 4.10 



where = /c^ + . In the first equality, we have substituted /c = and used the fact 
that the potential energy of thermal logarithmic potential is small compared to one of 
thermal mass potential, o?T'^ <^ T^. In the second equality, we have used 0c = T/X. This 
shows that no non-perturbative particle production occurs while the oscillates with the 
thermal mass for A -C a. 

Therefore, for A ^ a, the thermal potential derived in Sec. [3] is applicable at all times. 



(ii) A > a 

In this case, the oscillation frequency is much larger than the thermalization rate. Hence, 
the propagators of x cannot be considered as thermal ones and have to be treated as 
dynamical ones. There are two dominant effects that produce the x particles: production 
from thermal plasma and non-perturbative production. The essential difference from the 
case (i) is that even if the produced x particles become heavy due to the amplitude of 0, 
the produced particles may not decay and survive at each oscillation of 0, depending on 
the balance between the oscillation time scale of and the decay rate of x- The survived 
heavy particles significantly affect the motion of 0, since they form a linear potential 
effectively: 
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as in the case with the moduh trapping [80] . The motion of is determined by the comph- 
cated balance of parameters, and so, in the following, we will not consider the case (ii) for 
simplicity. The complete analysis of the trapping effect in the presence of thermal plasma 
will be performed elsewhere. 



5 Dynamics of oscillating scalar field 

Basic ingredients for analyzing the scalar dynamics have been presented above. Let us 
follow the evolution of the scalar field. First, we divide three cases depending on which 
term in the effective potential determines the oscillation epoch of : zero-temperature 
mass, thermal mass or thermal log. 

• Case (a) : begins to oscillate with thermal logarithmic potential if 



(j)^ < (f)^ = aTnJ — and A^f > TrM^^ (5.1) 

• Case (b) : begins to oscillate with thermal mass if 

• Case (c) : Otherwise, (f) begins to oscillate with zero-temperature mass. 

Note that we assume that x particles are absent initially in the ). Otherwise, 

would feel correction to the effective potential. This assumption might break down if 
heavy x particles are substantially produced by the direct decay of inflaton or the inflaton 
preheating process. Whether this occurs or not depends on the inflation model, and hence 
we simply assume the absence of x in the case (a). 

We take parameters so that T^/Mp < is satisfied in all the following analysis. 
Therefore, in any case, begins to oscillate during the inflaton oscillation dominated 
phase. The Hubble parameter at the beginning of oscillation, Hos, is then given by 

( a^TlMp/(j)j for case (a). 
Has ^ I {X^TlMpf'"^ for case (b), (5.3) 

for case (c). 

Fig. [T] shows contours of Hos/m^ on (0j,A)-plane for a = 0.05, Tr = lO^GeV and = 

1 TeV (top) and for = lO^TeV (bottom). Regions (a) - (c) correspond to the cases 
(a) - (c) described above. Note that it is possible that even if begins to oscillate with 
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Figure 1: (Top) Contours of Hos/m^ for a = 0.05, Tr 
(Bottom) Same as top panel, but for = 10^ TeV. 



10^ GeV and m^, = 1 TeV. 
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zero-temperature mass term, thermal effects become dominant thereafter. Conversely, the 
zero-temperature term eventually becomes dominant even if thermal effects are important 
at the beginning of oscillation. These facts make the scalar dynamics quite complicated. 

The evaluation of decay or evaporation epoch of is even more complicated due to 
the kinematical condition including the field-dependent and temperature-dependent mass 
of X well as the existence of thermal dissipation rate. In order to follow the dynamics 
of oscillating scalar field, it is convenient to study averaged quantities with a time interval 
which is longer than the oscillation period but shorter than the Hubble time scale. We 
summarize useful equations in Appendix El The averaged effective dissipation rate, F^*^, 
at an each averaged amplitude regime, which is defined in Appendix [B], is summarized 
as follows. 

• \i (j) oscillates with thermal log potential (^Ti^g ~ the dissipation is caused 
by scattering with gauge bosons in thermal bath. Notice that the (f) decay into gauge 
boson pair is kinematically forbidden. The dissipation rate is evaluated as 

Tf^ha-T'^^^. (5.4) 

We have shown that there are no efficient non-perturbative production in this regime 
for A < a. 

• If oscillates with a thermal mass term {ra^^ ~ AT), we need to distinguish two 
cases : A < a and A > a. In the former case, the oscillation is adiabatic with respect 
to thermal relaxation rate and the dissipation is caused by scattering with x particles 
in thermal bath. In the latter case, can kinematically decay into x pair. However, 
the dynamics is so so complicated and we do not consider this case. Thus we have 

Vf^X^aT for A0 < aT, (5.5) 



As a rough approximation, we simply extrapolate this result to aT < Xcf) <^ T|*^ 



We have shown that there are no efficient non-perturbative production in this regime. 

If oscillates with a zero-temperature mass term (m^^ ~ ^<i>)-, "we have several 
situations. First, consider the case of m^g -C aT. In this case, the dissipation 
coefficient can be evaluated similarly, since the motion of is adiabatic with respect 
to the thermalization time scale of hot plasma. 



\ha^T^/{^<P,^) for T<A0, 
A^aT for A0 < aT. 



"^^^ For numerical simplicity, wc do not use the dissipation coefficient computed from Eq. p.30p which 
is applicable at aT < A(/) ^ T. 
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As a rough approximation, we extrapolate these results to aT < A0 < T j^^l Then 



the threshold 0th can be estimated as gT/X if the non-perturbative production is 



absent or Ym00/A, where the adiabaticity is broken, if it is present. Next, consider 

the case of m^jj ^ gT. In this case, F^^ can be estimated with neglecting the finite 
density correction to the dispersion relation of x 
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leff 



[ a^m^^/ (00th) for < A0, 
IA^m(^/(87r) for X(p <^ Tiifj,. 



r; 



The threshold 0th is evaluated as ym^cp/X. Though we do not calculate at the 

regime aT < m^g < gT actually, we simply extrapolate these results as approximate 
ones. As pointed out in Sec. 14. the non-perturbative particle production gives the 
effective dissipation rate as 

A^m ~ f g'^T'^ 1 
^ - for A0 ^ max < m^, > (5.8) 



Given above formulae, we can trace the evolution of scalar field oscillation. One more 
complexity arises from the fact that the time dependence of quantities T, and so on 
changes before and after the reheating. After the reheating, can again dominate the 
Universe depending on parameters. The evolution equations of the system are given by 

+ (3if + r<^(0; T))0 + ^^^^^^ = 0, (5.9) 

pr+4:Hpr = TinlPinl + T^p^, (5.10) 



1 

'p 

.2 J,2 



H =^^ipmf + P^ + pr), (5.11) 



where V := m^0 /2 + V^g, pr = {n g^/30)T is the radiation energy density at the leading 
order in g and A, and pinf is the infiaton energy den sity and Finf = a/ {TT'^g^./90)T^/Mp is 



the infiaton decay rate. The energy density of q^lj is given by 



l02 + 1^2^2^ = + (5^12) 

In the second equality, we have used the virial theorem (See also Appendix [B]). Here 
denotes the time average with a time interval which is longer than the oscillation 



See footnote [#TT| 

As can be seen from Eqs. p.36|) and p.35p . there are Pauli blocking or Bose enhancement factors, 
though we will neglect them for simplicity in the following. 

The energy density should not be confused with the free energy density. 
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period but shorter than the Hubble time scale. Validity of the last term in Eq. (15.101) is 
restricted to the case when the oscillates with the zero temperature mass term, however, 
it is sufficient for practical use because the dominant change of temperature due to the 
loss of 0's energy occurs when the ^'s energy dominate the universe, that is, when the (f) 
oscillates with the zero temperature mass term. Following scaling relations will be useful 
(See Appendix [Bl) : 

for zero temperature mass, 

for thermal mass, (5.13) 
for thermal log. 

The scalar field is expected to evaporate when H = r^(T, 0) is satisfied. In order to 
see how the 0's amplitude evolves and to clarify thermal effects on the evolution of 0's 
amplitude, it is convenient to consider the quantity 




_ m20V2 

Krh = 



Pinf + Pr 



(5.14) 



In Fig. [2] we have plotted contours of for a = 0.05, Tr = 10^ GeV and = 1 TeV 
(top) and = 10'^ TeV (bottom). In the region labeled by "0-domination" , the energy 
density dominates the Universe before it decays. It is seen that the behavior is highly 
non-trivial. This kind of complicated structure would not appear without taking thermal 
effects into account. In order to see what is happening, let us see some typical cases in 
these parameter regions. 



5.1 Oscillation with thermal log 

After begins to oscillate with thermal logarithmic potential, the amplitude of 0, decreases 
as oc a~^/^ while the temperature decreases as T oc a~^^^. As described in Sec. 13.11 in 
this regime x has a large mass and is decoupled from thermal bath. By integrating out 
the heavy x field, has an effective interaction with gauge fields as £ ~ 6(j)FF/(j). Note 
that the effective mass of may be estimated as rrtg^ ~ aT"^ /(j) <^ gT, hence the decay 
into gauge bosons are kinematically forbidden. However, receives dissipative effects from 
thermal bath and the dissipation rate is given by F^''^ ~ ha^T^ / ip^. 

Fig. |3] shows time evolution of various quantities as a function of Hubble scale. At 
H = Hos ~ 6 X 10^ GeV, begins to oscillate with thermal logarithmic potential. As the 
amplitude decreases, the dissipation rate decreases F^^ oc H^^^ and it becomes equal to 
the Hubble parameter H at H^ec ~ 10^ GeV, where is expected to evaporate. In this 
case, therefore, the coherent oscillatio n so on disappears after the onset of oscillation due 
to the interaction with thermal plasma J^^^ 

'^^^ It may be the case that (p coherent oscillation deforms into non-topological solitons, oscillons, before 
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<f i [GeV] 

Figure 2: (Top) Contours of for a = 0.05, Tr = 10^ GeV and = 1 TeV. (Bottom) 
Same as top panel, but for = 10^ TeV. In the region labeled by "0-domination" , the (p 
energy density dominates the Universe before it decays. 
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Figure 3: Time evolution of various quantities as a function of Hubble scale for (A, (pi) = 
(2 X 10-3, IQis GeV). We have taken = 1 TeV and Tr = 10^ GeV. 



5.2 Oscillation with thermal mass 

Next, we consider a typical case where thermal mass plays an important role. Fig. H] 
shows time evolutions of various quantities as a function of Hubble scale for (A, (pi) = 
(10-^ 10^^ GeV), = 1 TeV and Tr = 10^ GeV. As can be seen in Fig. [H (p begins to 
oscillate with a thermal mass for this parameter choice at if = Hos ~ 10^ GeV. As the 
temperature decreases, thermal mass of (p also decreases as m^g = AT oc H^^^ in the 
matter dominated era and m^g oc H^^^ in the radiation dominated era after the reheating. 
At the temperature T ~ m^/X = 10^ GeV, or ~ lO^^GeV, the zero-temperature mass 
begins to dominate. Since <^ aT at this stage, the motion of (p is adiabatic and the 
main dissipation effect comes from the scattering with x in thermal bath. The dissipation 
coefficient is given by ~ A^aT. Thus (p evaporates at if = H^ec ~ 10~^ GeV, where 

iidcc ~ A^a^Mp. (5.15) 

The temperature at the evaporation is then estimated to be Tjcc ~ 10'' GeV, consistent 
with the assumption that it occurs after the reheating : Tjec < Tr. 

the evaporation. Since the time scale of the development of spatial instability is comparable to the 
evaporation time scale, we avoid definite conclusion about this issue. The case of oscillon formation will 
be discussed in the Appendix |^ 
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Figure 4: Same as Fig. [3l but for (A, 



;lO"^10^^ GeV). 



5.3 Oscillation with zero-temperature mass 

Let us consider the case where begins to oscillate with zero-temperature mass. In the 
top panel of Fig. |5l we show time evolutions of various quantities as a function of Hubble 
scale for (A, 0i) = (IQ-^, 10^^ GeV), = 1 TeV and Tr = 10*^ GeV. After the reheating, 
soon dominates the Universe aX H ^ 10~^ GeV. Since the amplitude is so large, x obtains 
a large mass and decouples from thermal bath. By noting that decay into gauge bosons 
is kinematically forbidden due to large thermal mass, the main dissipative effect comes 
from the scattering with gauge bosons in thermal bath. The dissipative coefficient is given 
by r,^ ~ ba^T^ / {(j)(j)th) ■ These thermal particles scatter off and as a result, evaporates 
at i7 ~ 10"^ GeV. 

In the middle panel of Fig. [5l we show time evolutions of various quantities as a function 
of Hubble scale for (A, 0^) = (10-^ 10^^ GeV), = 1 TeV and Tr = 10^ GeV. In this case, 
A is so small that the dissipation rate is suppressed compared with the above case. After 
the Hubble parameter decreases to if ~ 10~^ GeV, the amplitude of becomes small so 
that the decay into x P^-ir is accessible. Finally radiation generated from dominates 
over the inflaton decay products. The decay rate at this stage is given by F,^ ~ \'^m^/{8n). 
Therefore, decays at if = ifdec ~ A^m<^/(87r) ~ 10~^^ GeV. 

In the bottom panel of Fig.[5l we show time evolutions of various quantities as a function 
of Hubble scale for (A, 0^) = (10-^ 10^^ GeV), = 1 TeV and Tr = 10^ GeV. The 
situation is similar to the former case, but in this case the amplitude of is so small and 
never dominates the Universe. Thus decays at if = ifdec ~ A^m(^/(87r) ~ 10~^^ GeV, 
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when the cosmic temperature is around T ~ 100 GeV. 

6 Conclusions and discussion 

In this paper we have extensively studied the scalar dynamics in the early Universe taking 
account of the effects of thermal environment. Despite the fact that the early Universe 
after inflation may be filled with high-temperature thermal plasma, thorough analyses 
on the scalar dynamics in such an environment were missing in the literature. Needless 
to say, scalar fields play important roles in cosmology. Infiaton explains the primordial 
inflation and the primordial density perturbation, and curvaton may also be responsible 
for the generation of density perturbation. Affleck-Dine fields in the MSSM can create 
the baryon asymmetry of the Universe. SUSY breaking fields may have significant effects 
on cosmology. Similar long-lived scalars often appear in extensions of the SM. All these 
scalar fields generally have large energy densities so that they must decay or evaporate at 
an appropriate epoch in order not to disturb the success of standard cosmology. Therefore, 
they necessarily couple to SM particles directly or indirectly through some intermediate 
states. 

Based on these observations, we considered a model in which a scalar (p couples to 
fermions x through a Yukawa coupling, which then interact with thermal bath. This 
simple model captures essential features of realistic models. We have consistently taken 
into account various effects : thermal modification on the effective potential of 0, thermal 
dissipation of </>, non-perturbative particle production and the formation of non-topological 
solitons. It is found that even in this simple class of models, the scalar dynamics is 
so complicated that a naive estimation neglecting these effects is not allowed in broad 
parameter spaces. In particular, it should be noticed that thermal dissipation, which can 
be interpreted as the evaporation due to scatterings with particles in an environment, 
often plays a dominant role in determining the dissipation epoch of 0. 

Finally, let us mention what we have not included in our analyses. In this paper we 
have followed the evolution of only the zero-mode of coherent oscillation. However, it 
is expected that fluctuations around the zero-mode, or the particle-like excitations of 
can be important if thermal effects play a dominant role for the evaporation of 0. This 
can in principle be traced by following the evolution of the two-point functions of 0. We 
have also restricted ourselves to the case of A < a, because otherwise the whole dynamics 
including x particle production is complicated. We will return to these issues in some 
concrete models elsewhere |91j . 
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A Formation of non-topological solitons 

We have seen that a scalar field (f) may begin to oscillate in the thermal logarithmic 
potential, depending on parameters. It is known that, however, a coherently oscillating 
scalar field with a potential flatter than the quadratic potential exhibits an instability 
that results in formations of classical lumps. These are called oscillons and studied in 
the context of inflaton fragmentation (HH [Ml |83l [H [851 186] . In Ref. [87], it is dubbed 
as "I-ball" in analogy with the Q-ball whose stability is ensured by a conserved U(l) 
charge. Similarly, the stability of I-ball is guaranteed by the existence adiabatic invariant 
in association with the scalar dynamics, which is often denoted by /. Once the scalar 
field fragments into oscillons, their subsequent cosmological evolutions may be modified. 
Although it may have non-trivial consequences, the dynamics is hig hly non- linear and it is 



difficult to derive robust conclusions without numerical simulations^ Thus here we will 
only shortly see possible effects of oscillon formation in this Appendix. 

It is noticed that even for the zero-temperature mass dominant case, the oscillon for- 
mation may take place because the mass receives radiative corrections from the fields that 
couple to 0, and it may modify the mass like m'^{(j)) ~ '^'^{1 + eln(</))), where e represents 
a one loop factor including coupling constants. If e < 0, the oscillon formation may take 
place. We do not specify here the model-dependent radiative correction and simply con- 
sider the finite-temperature mass dominant case in which the oscillon formation cannot 
be avoided. 



A.l Oscillons 

First we summarize the properties of the oscillon. The most discussion follows from that 
performed in the context of Q-balls [261 [23 [281 [291 [301 [3l] • 

As the (j) begins to oscillate with thermal log term, low frequency modes exhibit insta- 
bilities to form the oscillon. The most important mode is ~ Hos ~ « ^TlMp/(j)l and the 
corresponding oscillon solution has a typical radius of i? ~ 1/Hos- According to Ref. [87], 
the oscillon configuration corresponds to the bounce solution with fixed adiabatic invari- 
ant, /. In this case, it corresponds to the particle number density : Trf^cjP'. Therefore, the 
total charge of an oscillon is estimated as 

I - I^^^H^A ^ ^\ — ] , (A.l) 



OS 



#16 xj]-|2g scale required for the development of oscillon configuration is more or less close to the dissipation 
time scale in the case of thermal logarithmic potential. This makes it difficult to estimate the probability 
for the oscillon formation. 
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<fi [GeV] 

Figure 6: Oscillons may be formed in regions inside the line for = 1 TeV (black thick) 
and = 10^ TeV (blue thin). We have taken Tr = 10^ GeV. 



where /3 is a numerical factor which represents the delay of oscillon formation from the 
epoch of H = Hos- 

Notice that even if (j) begins to oscillate with zero-temperature mass term, thermal 
correction may come to dominate the effective potential depending on parameters. In this 
case the "delayed" -type oscillons may be formed [32]. Since the time scale of oscillation 
is much faster than the Hubble scale and the dissipation rate, the oscillon formation may 
be inevitable in this case. 

Fig. [6] shows the parameter regions where oscillons may be formed for = 1 TeV 
(black thick) and = 10^ TeV (blue thin). We have taken Tr = 10^ GeV. Inside these 
lines, (f) experiences the era of thermal log dominated regime before it evaporates and 
hence oscillon formation potentially affects the fate of 0. 

A. 2 Decay of oscillons 

At the classical level, oscillons are regarded as stable objects [88l[89], although they can 
decay at the quantum level [90]. In the present situation, the problem is not so simple 
because the oscillon solution ceases to exist as the cosmic temperature decreases and the 
thermal potential becomes subdominant. Thus oscillons are not absolutely stable in our 
setup. Once oscillons are formed, the final decay temperature of (p condensation may not 
be estimated by using the perturbative decay rate of 0. There are a couple of effects that 
makes oscillons unstable, (a) Perturbative decay of (p. (b) Deformation of the oscillon due 
to the Hubble expansion. 
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A. 2.1 Perturbative decay 



Inside oscillons, (f) has a large amplitude of (pi and x obtains a large mass. As noted before, 
the preheating is not efficient if (f) oscillates around the thermal logarithmic potential. 
Therefore, we focus on perturbative decay of (f) inside and at the surface of the oscillon. 

Denoting the scalar field configuration of an oscillon by (f){r) where r is the radius 
measured from the center of oscillon, the decay rate depends on r. For r < Tc where is 
defined as (/>(rc) = 0c = T/\, the x fi^ld decouples from thermal bath. For r > Vc, x is so 
light it can take part in thermal bath. Thus the decay rate is roughly expressed as 

X^aT for r > Vc- 

This shows that the decay of oscillon occurs most efficiently from the surface at r ~ r^. 
Hereafter, we make use of the approximation of ~ R, since the typical length scale of 
the change of the field value is R. The oscillon charge evaporation rate is then roughly 
approximated by 

- f = / rfr4vrrHV(r)^r,(r) ~ ^Xa^T^ j . (A.3) 

The effective oscillon decay rate is given by 

r--l^ o^AoT^ i/3-^aXH{^/(j),) for T > Tr 

' ~ Idi^ ~M^Jc ^ \/3-iaAff(T/TR)2(0/0,) for T < Tr. ^ ' ' 

Thus this implies that the Tj is smaller than H unless /3 is so small or (0/0c) is so large. As 
the temperature decreases, the zero-temperature decay rate becomes dominant and then 
it may exceed the Hubble rate. Before then, the oscillon may be deformed as explained 
below. 



A. 2. 2 Deformation 

Since the cosmic temperature decreases as the Universe expands, the scalar potential will 
finally be dominated by the zero-temperature one. Then a stable oscillon solution does 
not exist and oscillons are expected to disappear. This happens when the temperature 
drops down to T ~ Tdcf = {rn^cpi/ aY^"^ . In terms of Hubble parameter, it corresponds to 
H ~ ifdcf given by 



E^,, = { a^TjMp /7os . (A.5) 



— — for Tdcf < Tr 

aMp 
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After the oscillon deformation, dissipates its energy into thermal bath in the way de- 
scribed aheady. If the oscillon survives long enough, it is possible that oscillons dominate 
the Universe since the oscillons behave as matter. 

B Evolution of oscillating scalar fields 

In order to follow the evolution of oscillating scalar fields, it is convenient to consider 
averaged quantities with a time interval which is longer than the oscillation period but 
shorter than the Hubble time scale. In particular, we will derive the formulas which 
describe the evolution of 0's amplitude and the evolution of 0's energy density, and clarify 
the effects of dissipation on them. 

B.l Relevant equations 

Let us consider scalar fields oscillating in the potential V{(f); T) = m?^(j? + Kff(0; T), which 
includes thermal effects. In the presence of dissipation. The equation of motion is given 
by 

0+(3i? + r<^)0 + |^ = O. (B.l) 

From this, we immediately find the following relation between the kinetic energy K = 0^/2 
and the potential energy V : 

|(if+V') = -(6i/ + 2r,)A- + I(Tf), (B,2) 

This is an exact relation. This means that + is changing slowly 0{H,r^), and 
approximately conserves with time. Therefore, we can take time average of this equation 

as 

j^(KTV) = -QH(K) - 2(TV^) + ^ (t^^^ , (B.3) 

with a time interval which is longer than the oscillation period but shorter than the Hubble 
time scale. In this section, the time average is denoted as (^^). 

First, we make use of the virial theorem for deriving the averaged motion of (p. By 
noting that 2K = d/dt{(f)<p) — (jxp, we have 

(27f> = (*§^). (B.4) 

with neglecting 0{H^,r^H) terms. This is the virial theorem for an oscillating scalar 
field. 
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Next, let us evaluate {2T^K) with some examples. 

(i) Zero temperature mass at the large amplitude regime: In this case, the dependence 
of dissipation coefficient is given by 1/0^ above the threshold value 0th which should be 
determined case by case. Then one finds 



ml \ , f , 1 1 



(r^</) ) oc ( ^^^2^^^ j "^'^ Jt^^'^^ tan^ m^r "^'^tth 

~ -^(5^) (B.5) 

0th0 

Here we have used (p oc (j) sin m^t, (pth ~ (pm^tth and the virial theorem {K) = {V). 

(ii) Thermal logarithmic potential: In this case, the dissipative coefficient is given by 

T^c^ba^ — . (B.6) 



Hence we obtain 



(r<^02) = ^^2^3 (^r^^ ^^2^3 l ^^^^-^y ^ (02). (B.7) 

Here we have used the energy conservation 0^/2 = a^T^ ln(0^/0^) and the virial theorem 
(02) = 2a^T^. Therefore, the averaged dissipative coefficient is obtained 



/ ln(0?/0' 



^^2y3/ ln(0|/0^) y ^3_g) 



In the following, we denote the averaged dissipative coefficient as F^^ collectively. Thus, 
Eq. ( ]B.3|) can be expressed as 



d 



{K+V) = -{QH + 2Tf){K) 



dV 



T- 



(B.9) 



From Eqs. f lB.Op and (lB.4p . we obtain a time-averaged relation for the scalar potential as 



dV 



(6i/ + 2F-) 0— +- T— . 



2T IdV 



T 



dT 



(B.IO) 
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B.2 Evolution of scalar amplitudes in general effective potential 

Let us derive the evolution of scalar amplitude with some examples using Eq. ( IB.IOI) . 

(i) Zero temperature mass : V = {l/2)rn?^(f? . Substituting it into ( IB. 101) . we obtain 

j^m = -{^H+Tf)m. (B.ii) 



Writing (0^) = 0^, we have (j)^ (x a ^ aX the regime F^^ <^ H 



#17 



ii) Thermal mass : V = (1/2)A2tV • Substituting it into flR10|) . we obtain 



|(0^) = -(^3if + P/ + ^j(02). (B.12) 

Thus we have 0^ oc a~^T~^ at the regime <^ H. 

(iii) Thermal log : V = a'^T^ln{\'^(j)^ /T^). Substituting it into f lRTOj) . we obtain 

— (T4(l + ln(/.VT2)) = -(6i7 + 2rf )T^ + — (T4 1n02/T2) T\ (B.13) 

at V rp 

The l.h.s. can be written as 

— {T\l + In 02)) + T'-U^ - —T' (B.14) 



where we have defined (In 02) = In^^. Thus we have 



6i7 + 2r/ + — U^. (B.15) 



Thus we have 02 (x ^j^g regime F^''^ ^ H. 

B.3 Evolution of the energy density of scalar field 

Let us derive the evolution of 0's energy density, which is defined as 



P4> = (^0^ + ^"^10^) = (K) + ^mj(02). (B.16) 



The definition of (j) here is a bit different from that used in SecH) The difference, however, is at most 
0(1) and we do not distinguish them since our discussions do not require 0(1) accuracy. 
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To make our discussion concrete, we will study the evolution of energy density with some 
examples in the following. Note that the evolution of the latter term, i.e. the potential 
energy, have been obtained in the last section. 



(i) Zero temperature mass : V = (l/2)m^0^. Using the virial theorem: Eq. ( IB. 41) . we 
obtain 

(K) = \ml{W). (B.17) 
Thus, the evolution of {K) is the same as one of (0^) given by Eq. ( IB. Ill) : 



j^P, = j^{2K) = mj|(02) = _(3if + Tf)p,. (B.18) 



(ii) Thermal mass : V = (1/2)A^T^0^. The virial theorem implies 

(K) = Ix'T^-^). (B.19) 

As can be seen from this equation, if the oscillates with the thermal mass potential, 
then {K) ^ m^^{(fP')/2 is satisfied. Differentiating this with respect to time, we find 

j^{K) = -(^H^Tf-p^{K). (B.20) 
And we also have Eq. ( ]B.12p : 



I = - 3if + rf + ^ /im|02 ) . (B.21) 



(iii) Thermal log : V = a^T'^ln(A^(/)^/T^). The virial theorem implies that the kinetic 
energy is independent of the dissipative coefficient: 

(K) = a^T^. (B.22) 



Hence we have 



and Eq. flB.lSP means 



I (^^0<^7 = -[qH + 2Tf + ^ j . (B.24) 

Note that the kinetic energy is larger than the potential energy if the oscillates with the 
thermal logarithmic potential: {K) ^ (m^02/2). 
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C Closed Time Path formalism 



In this section, let us briefly introduce the basic ingredients of Closed Time Path formalism. 
Though equations and formulas shown in the following can be found in [601 l6T| [621 [7] and 
references therein, we will summarize them for the sake of readers. 



C.l Closed Time Path formalism 

To follow quantum dynamics, we are often interested in the time evolution of expectation 
value of operators for a system described by a density matrix p at an initial time. Such 
expectation values can be calculated in terms of so-called Closed Time Path (CTP) or 
in-in formalism (See [601 161] and references therein). The evolution of expectation value 
of a Heisenberg operator can be written in terms of CTP as 



{Onit)) = tr 



pTcexp^^-z j dt'Hiit')^ Oiit) 



(C.l) 



where the H and I subscript denote the Heisenberg and Interaction picture respectively, 
the time integration is performed on the CTP contour C, and Tc denotes the contour C 
ordering. The Schwinger-Keldysh propagator is defined as a connected two point correlator 
with the CTP contour: 

G(x, y) := (Tc (^(x)(^(y)) - 0(a;)0(y) (C.2) 

S{x,y):= {Tci^{x)i,{y)), (C.3) 

where 0(a;) := {(p{x)). Two convenient combinations of propagator are introduced: the 
Jordan propagator (or the spectral function) and the Hadamard propagator (or the sta- 
tistical function). 

Gj{x, y) := mx), m]); Gh{x, y) ■= ({^(x), ^(y)}) - 2<P{x)<P{y) (C.4) 
Sj{x,y) := {{^{x),$iy)}); 5^,(x,y) := ([^'(x), ^(y)]). (C.5) 
Using these propagators, one can express the Schwinger-Keldysh propagator as 

Gix, y) = ^ [Gh{x, y) + sgnc-(xo, yo) Gj(x, y)] , (C.6) 

S{,x,y) = ]^[SH{x,y) + sgncixQ.yo) Sj{x,y)] (C.7) 

where the sign function sgn^; is defined on the contour C. It is often convenient to define 
the retarded and advanced propagators: 

G'rct/adv(a;,y) := ±i9{±XQTyo)Gj{x,y) (C.8) 
5'rct/adv(a;,y) := ±i9{±XQTyo) S.j{x,y). (C.9) 

Note that if a system has a spacial translational invariance, all the above two point 
correlators only depend on the difference of spacial coordinate, x — y. 
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C.2 Thermal Equilibrium 

In general, the Hadamard and Jordan propagators are independent. However, if some fields 
are in thermal equilibrium (See [62117] and references therein), their Hadamard and Jordan 
propagators are related through the Kubo-Martin-Schwinger (KMS) relation [71] . Due to 
the translational invariance of thermal equilibrium system, all the two point correlators in 
thermal equilibrium depend only relative coordinates and it is convenient to consider the 
Fourier transformations of them: 

(x - = / ^ (A;). (C.IO) 

The KMS relation implies 

G'^iu, k) = [1 + 2Muj)] p^(a;, k), (C.ll) 

S'^iu, k) = [1 - 2Muj)] p^iu, k), (C.12) 

where /b/f is the Bose-Einstein/Fermi-Dirac distribution and p^^^ denotes the spectral 
densities, defined as p^/^{u,p) := G^j / S^J'{uj,p) respectively. 

For simplicity, the "th" superscript is suppressed in the following. Eqs. (IC.Sp and (1C.9P 
imply the following relations: 

G'rct/adv(PO,P)l ^ pY /" ^^0 ^ |Gj(A;o,P)| ^ i_ |G'j(PO,P)| 

'S'rct /adv (Po,p) J J 27r j>o - A;o 1 5j(/co, p) J \ Sj{po,p) ] ' 

The retarded and advanced self energies are defined as 

Sret/adv(a;) J [Sj(x)J 

and these imply the following relations: 

nrct/adv(Po,p)l ^ pY /" ^fep 1 |nj(/i;o,p)| ^ i_ |Hj(po,p)| 
-I ret /adv (Po,p)I J 27r po - ^0 I Sj(A:o,p) J 2[Sj(po,p)l 

If a real scalar field (f interacts with thermal bath by Ci^t = fO, then the self energy 
of (f is given by 

Hj(x) = ([6(x),6(0)]) (C.16) 

at the leading order in Ci^t- Here we have assumed that the back reaction on the thermal 
bath is negligible. In this case, the real and imaginary parts of retarded self energy can 
be expressed as 

3?Hret(po,p) = PV / — (C.17) 

c-TT ^ ^ nj(po,p) 

3nret(po,P)= ^ • (C.18) 
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In an interacting theory, these spectral densities can be expressed as 

1 1 



P^(Po,p) 
P^{Po, P) 



_m2 - p2 + riretbo, P) 771^ - + Ti^dv{Po, P) . 



(C.19) 
(C.20) 



where 11 and S are the self energies for boson and fermion, and m/5 denotes the Dirac 
mass (See [73] for Majorana fermion). 

C.2.1 Breit-Wigner Approximation 

If the spectral densities well concentrate around the poles, we may approximate them by 
the Breit-Wigner form. For a real scalar field, the Breit-Wigner form is given by 



P (Po,p) 



2porp 



12' 



(C.2I; 



m2 + p2 JJIIretl^^p, p) and Tp = -55ni.ct(^^p, p)/^^p- 



where Qp = 

As shown in Ref. |6^, the Breit-Wigner form of spectral density for fermion with the 
vanishing mass term can be expressed as 



P^(Po,p) 



(70 - q • 7) 



2 [[po-%? + r^py^"" " " [Po + %? + r^py^ 

where the real and imaginary parts of the pole are determined by 



(70 + q ■ 7) 



(C.22) 



±p-AS{nt){niTp)-B^{ni) 

and the wave functional renormalization is given by 

d 



(C.23) 
(C.24) 



1 



dco 



UJ 



(C.25) 



From Eq. (1C.15I) . the self energies can be decomposed into Srot/adv = S =1= iS j/2, and can 
be expressed as 



E(po,p)= -A^{po)^-B^{po)-fo 
Sj(po, p)/2 = A'pipo) i> + fi^(po) 70- 



(C.26) 
(C.27) 
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On the other hand, as shown in Ref. [TO], if the zero temperature mass term is much 
larger than the typical thermal mass scale gT, then the Breit-Wigner form of spectral 
density for fermion can be approximated by 

where the real and imaginary parts of the pole are determined by 



ml + p2 (C.29) 

rp/2 = — [ml + C^K)) + B^Mu^] . (C.30) 

Up 

In the presence of Dirac mass term which breaks the chiral symmetry, the self energies 
can be expressed as 

S(po, P) = - A^ipo) i> - B^ipo) 70 - C^iPo) mn (C.31) 

Sj(po,p)/2 = Alipo)^ + Bl,{po)^o + Clipo)mD. (C.32) 
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